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Abstract. An update is presented to the software package SeBa (Portegies Zwart and Verbunt 
[1], Nelemans et al. [2]) for simulating single star and binary evolution in which new stellar 
evolution tracks (Hurley et al. [3]) have been implemented. SeBa is applied to study the population 
of close double white dwarf and the delay time distribution of double white dwarf mergers that may 
lead to Supernovae Type la. 
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SEBA - A FAST STELLAR & BINARY EVOLUTION CODE 

We present an update to the software package SeBa (Portegies Zwart and Verbunt 
[1], Nelemans et al. [2]) for simulating single star and binary evolution from the ZAMS 
until remnant formation including processes as mass transfer phases, common-envelope 
phases, magnetic braking and gravitational radiation. Previously stellar evolution has 
been based on evolutionary tracks described by analytic formulae given by Eggleton 
et al. [4] (hereafter EFT), but in the new version it is based on Hurley et al. [3] (hereafter 
HPT). In this research SeBa is used to study close double white dwarfs. A comparison 
between simulations using the EFT and HPT stellar evolution tracks shows that the over- 
all double white dwarf populations are similar, see Fig. 1 . In the models presented here 
we assume solar metalicity and a 50% binary fraction. Inital masses of single stars and 
binary primaries are distributed according to the Kroupa ([5]) IMF. Secondary masses 
are drawn from a flat mass ratio distribution. Initial orbital parameters are distributed ac- 
cording to a thermal eccentricity distribution and oc l/a (Abt [6]). Differences between 
the populations are traced back to differences in the evolutionary tracks that affect the 
stability of mass transfer (at M1/M2 ^1.5) and the size of core masses of AGB stars (at 
M1/M2 < 1). 
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FIGURE 1. Simulated distribution of population of visible double white dwarfs as a function of orbital 
period and mass ratio, where mass ratio is defined as the mass of the brighter white dwarf devided by 
that of the dimmer white dwarf. Left the stellar evolution tracks according to EFT are used, right HPT. 
The intensity of the grey scale corresponds to the density of objects. Observed binary white dwarfs are 
overplotted with filled circles. 



DELAY TIMES OF DOUBLE WHITE DWARF MERGERS 

Close double white dwarfs tend to lose considerable amounts by the emission of gravita- 
tional radiation, which is not only interesting for future gravitational wave observatories 
as LISA (o.a. Evans et al. [7]), but also affects the binary system by decreasing the orbital 
period and possibly leading to a merger. The merger of two carbon-oxygen white dwarfs 
with a combined mass exceeding the Chandrasekhar mass limit is one of the possible 
scenarios for Type la Supernovae (SNe la). SNe la are very succesfuUy used as stan- 
dard candles on cosmological distance scales, but so far the nature of the progenitor(s) 
is unclear. 

The delay time distributions of these mergers for the EFT and the HPT stellar evolu- 
tion tracks for solar metallicities are very similar. The delay time distribution for both 
types of tracks show that these mergers are expected to take place in young as well as 
old populations. The time-integrated number of SNe la per unit formed stellar mass is 
4.3-10"'^ Mq-^ and 4.4- 10""^ Mq^^ when using the EFT resp. HPT tracks. From cluster 
SN la measurements, Maoz et al. [8] infer a value of 5.9- 10"-^ M©"^ for their 'optimal' 
model. The current merger rate with a Galactic star formation rate of Boissier and Prant- 
zos [9] is 1.4- 10^ yr^ for EFT as well as for HPT, where as the empirical Galactic SN 
la rate from Cappellaro et al. [10] is (4±2)-10~^yr~^ 
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FIGURE 2. Merger rate of double carbon-oxygen white dwarfs with total mass above the Chan- 
drasekhar mass as a function of delay time. Rates are in (100 yr)~^ per 10^^ MQofthe parent galaxy 
(SNuM). a) Black represents a simulation using the stellar evolution tracks according to HPT, grey EFT. 
b) Delay times using the HPT tracks are shown for three different prescriptions of the CE phase. In 
light- grey model aa, in grey model yy and in black model ya. 



COMMON ENVELOPE EVOLUTION 

Close double white dwarfs are believed to encounter at least two phases of common 
envelope (CE) evolution. In spite of the importance of the CE phase, it remains poorly 
understood. Several prescriptions for CE evolution have been proposed. The a formal- 
ism (Webbink [11]) is based on the conservation of orbital energy The a-parameter 
describes the efficiency of which orbital energy is consumed to unbind the common 
envelope according to: 

where Eorb the orbital energy and Egr is the binding energy between the envelope mass 
Ml and the mass of the primary M\. Egr is often approximated by: 



Egr= ^ T:'"" (2) 



GMiMi^env 

XR 

where R is the radius of the primary and X depends on the structure of the primary star. 
We assume aX = 2. 

In order to explain to the observed distribution of double white dwarfs Nelemans et al. 
[12] proposed an alternative formalism. According to this / formalism, mass transfer is 
unstable and non-conservative. The mass loss reduces the angular momentum of the 
system in a linear way according to: 

Jinit-Jfinal AM 

= y ,^ , ,^ (3) 



Jinit M1+M2 

where Jinu resp. J final is the angular momentum of the pre- respectively post-mass 
transfer binary and M2 is the mass of the secondary. We assume a value of 1.75 for 

7- 



In Fig. 2 we compare the delay time distribution for three different models of CE 
evolution. In model (xa the oc formalism is used to determine the outcome of every CE, 
in model 77 the /formalism is used every time. The preferred model ya is a combination 
of the a and 7 formalisms as in Nelemans et al. [12]. In systems with a giant primary 
and remnant or giant companion the CE phase is according to a, for a giant primary and 
a non-remnant or non-giant 7 is used. The 70c model shows high rates at all delay times, 
where as the rate for model aa decreases significantly at long delay times and for 77 
saturates at a level below 0.5 (100 yr)^ per 10^*^ MQat short delay times. As a result the 
current Galactic merger rates of model aa (5.9-10^^ yr^^) and 77(1.2-10^-^ yr^^) are 
~ 40% and ~90% respectively of that of model ya (1.4-10^ yr^)- The time-integrated 
number of SNe la per unit formed stellar mass is 4.4-10-'^ M© 1, 2.8-10-4 M© 1 and 
S.S IO-"^ M©-^ for model 7a, aa and 77 respectively. 
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